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ABSTRACT 

We present temporal and spectral properties of a unique X-ray dip in GX 301-2 as seen with 
Rossi X-ray Timing Explorer in May 2010. The X-ray pulsation from the source gradually de- 
clined prior to the dip, disappears for one spin cycle during the dip and is abruptly restored in 
the spin cycle immediately after the dip. Moreover, the phase-integrated spectrum of the source 
becomes softer before and during the dip and it quickly hardens again following the dip. Our 
findings indicate the fact that the mechanism for pulsations gradually turned off briefly and un- 
derlying dim and softer emission likely from the accretion column became observable in the brief 
absence of high level emission due to wind accretion. 

Key words. X-rays: stars - stars: pulsars: individual: GX 301-2 



1. Introduction 

GX 301-2 (also known as 4U 1123-62) is a High Mass X-ray Binary (HMXB) consisting of a 
neutron star and the early type companion, Wray 977 (BP Cru). The two objects are in an eccentric 
orbit (e = 0.462) with an orbital period of 41.5 days ( |Koh et aL] |1997| l. The spin period of the 
neutron star is 686 s (Sato et al. 1986[ l. The companion star, Wray 977 is often classified as a 



B1.2Ia star ( |Koh et al!l|1997[ ), while |Kaper et al.| ( ri995| l reclassified it as Bl lae-i- hypergiant at a 
distance of 5.3 kpc. The latter classification makes Wray 977 with a luminosity of 1.3 x 10^ and 
a mass of at least 48 one of the most luminous and most massive stars in our galaxy. The very 
high mass loss rate of M ~ lO^^Moyr"' gives rise to a very dense but slow (vco = 400 kms"') 
stellar wind. 

The neutron star in GX 301-2 accretes directly from the dense stellar wind, resulting in a highly 
variable X-ray light curve that is typical for wind accretion. X-ray luminosity changes by a factor 



of five within a time frame of an hour is common in GX 301-2 (Rothschild & Soong 1987 1. The 



orbital light curve shows a distinct profile: shortly before periastron, the source exhibits X-ray flar- 



ing activity, (the X-ray luminosity increases by a factor of ~25 Pravdo et al. 1995| . Around orbital 
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phase 0.2, the X-ray intensity reaches a minimum as the source passes through the dense inner 



regions of the stellar wind until it is almost eclipsed by the stellar companion (see e.g. Pravdo & 



|Ghosh| [200 1 1 |Leahy 1 12002| and references therein). Following the minimum, the X-ray luminosity 
increases again with a weak secondary maximum near apastron. This orbital modulation is ex- 
plained by [Leahy] ( |2002| l by the presence of a spirally formed gas stream flowing from the optical 
companion to the neutron star. The neutron star intercepts the gas stream twice per orbit: shortly 
before periastron where the neutron star runs along the path of the gas stream resulting in intense 
X-ray flaring activity and also near apastrom giving rise to the secondary peak in the orbital light 
curve. 

Earlier spectral studies of GX 301-2 could describe its X-ray spectrum with a power law mod- 



ified by an exponential cutoff at higher energies (White et al. 1983 i. Recent studies, however, were 
unable to model the spectrum of using simple models, instead they had to employ partial cov- 



ering models as the spectrum of GX 301-2 is heavily absorbed (see e.g. La Barbera et al. 2005 
Kreykenbohm et al.[[2004| with A^h ranging from lO^"* up to 2 x 10^"* cm"^ ([La Barbera et al.[[2005]) 



Around 6.4 keV a very strong iron fluorescence line with an equivalent width from 200 eV up to 
1800 eV is present (Leahy et al. 1988| l. At energies above 20keV, the spectrum is further modified 
by the presence of a cyclotron resonant scattering feature around ~ 37keV (Mihara 1995 1, which 



was shown to be strongly pulse phase dependent in energy and depth (Kreykenbohm et al. 2004 1. 

In this paper we uncover an X-ray dipping behavior that has seen in GX 301-2 for the first time. 
In the next section, we introduce the 2010 Rossi X-ray Timing Explorer (RXTE) observations 
of the source. In §3, we present our detailed temporal and spectral investigations of the RXTE 
observations on 2010 May 28 during which the source exhibited the X-ray dip and finally we 
discuss our results and possible mechanisms that could result in the observed behavior in §4. 



2. Observations and Data Processing 

We employ publicly available Rossi X-ray Timing Explorer (RXTE) observations of 
GX 301-2 which were performed in 17 pointings between 2010 May 26 and July 22 for a total 
exposure time of 184.7 ks (Program IDs: 95354 and 95428). Individual pointings had exposures 
ranging between about 1 and 34 ks. Those observations performed under the program ID 95354 
lasted over 10 ks each, providing a set of uniquely long and uninterrupted (by the satellite orbit 
around Earth) X-ray data of GX 301-2. 

In our study, we used the data collected with the Proportional Counter Array (PCA, Jahoda et 
al. 2006) only. Typically, one or two proportional counter unit (PCU) out of five available PCUs 
were operating during the observations. The combination of operational PCUs varied slightly be- 
tween observations. We employed all available data for the selected observations for timing studies 
while we used data collected with PCU2 only for time resolved spectral investigations to avoid 
calibration uncertainties (see §3.2). We processed RXTE data using tools of HEASOFT version 
6.9 and performed spectral analysis with XSPEC 12.5. In (Arnaud 1996). 



3. Data Analysis and Results 

We uncover an X-ray dip in the light curve of GX 301-2 as seen on 2010 May 28 for the first 
time. In Figure[T] we present the 3-25 keV light curve of the entire exposure span of about 12.4 ks 
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Fig.l. Light curve of RXTE/PCA observations of GX 301-2 on 2010 May 28 (ObsID: 95354- 
03-03-00) in the 3-25 keV range. The vertical dotted lines indicate the time interval of detailed 
temporal and spectral investigations presented in §3. 



RXTE pointing of the source. X-ray pulses with a period of 686 s are clearly evident. We find that 
the X-ray intensity of the source is smoothly decreasing in the episode leading to the dip which 
lasts about 1000 s (from MJD 55344.044 to 55344.055). X-ray intensity of the source gradually 
increases after the dip and exhibits typical fluctuations as seen from GX 301-2 . Our search in the 
light curves of all 17 RXTE pointings resulted in no additional dipping episode of the source. 

We present below the results of our detailed temporal and spectral investigations of RXTE 
observations on 2010 May 28 (corresponding to the orbital phase of ~0.7 using the ephemeris by 
Doroshenko et al. [2010]) to reveal the nature of the dipping behavior 



3.1. Temporal Properties 

In Figure|2j we show the 3-25 keV band light curves of 12 spin cycles of GX 301-2 in the episode 
leading to, during, and following the dip. Note that the vertical axis is plotted in logarithmic scale 
to better identify emission structures with low count rates. The pulse profile of GX 301-2 is typi- 
cally characterized by two broad structures; a main pulse and a slightly weaker second peak. We 
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Fig. 2. Spin phase light curves of GX 301-2 in the 3-25 keV band. Each individual plot is one spin 
cycle of the pulsar in a continuous time sequence within the time interval indicated in Figure [T| 
Notice that Y axes are in logarithmic scale. 

find that in the epoch leading to the dip, while the overall intensity of the X-ray emission de- 
creases, that of the second peak drops faster. Surprisingly, in the dip (cycle number 6), a weak 
emission structure at the spin phase of the second peak remains while the main pulse structure 
from GX 301-2 completely disappears. Following the dip, both emission structures are restored 
abruptly with the second peak being greater than or equally strong as the first one for four cycles 
(number 7-10). After that, the pulse shape of the source resembles its typical long term appearance 
again. 

We further investigated evolution of the pulse profile around the dip by means of disentangling 
the changes in the main pulse and second peak of the pulse profile. For each of the 12 cycles (Figure 
|2]) we fitted the pulse profile with a model consisting of a constant and two Gaussian components 
to account for the two peaks. As the constant level was not well constrained as it traded off with 
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Fig. 3. Evolution of integrated Gaussian functions 1 and 2 that are used to fit the main pulse and 
the second peak of the pulse profile. 

the Gaussian normalizations, it was fixed at the rate of 8.4 counts/s/PCU, which is the DC level 
in the dip. In order to ensure continuity between phase 1 and phase 0, which correspond to the 
same point, we fitted three versions of each Gaussian shifted by one period ahead and behind. In 
other words, each pulse was modeled with three Gaussians with a fixed spacing of one period, the 
same width and normalization. The fits are homogeneous across the cycles. Note however that the 
Gaussian fits to some of the cycles are not perfect (x^ ranging between 8.4 and 49 for 10 degrees of 
freedom), which is not surprising since the pulse profile of GX 301-2 is more complex than only 
two Gaussians can describe. Nevertheless, our main purpose is to evaluate time variation of the two 
broad pulse structures that can be tracked well with two Gaussians. We then computed the integral 
of each Gaussian (1 for the main pulse and 2 for the second peak) in the phase range [0,1] which 
we show in Figure|3]as a function of pulse cycle. It is important to note that both Gaussian integrals 
gradually decline during the cycles leading to the dip and similarly restored gradually following it. 

3.2. Time Resolved Spectral Properties 

We performed a time resolved spectral analysis in order to uncover the nature of X-ray emission 
from GX 301-2 prior to, during, and after the dip. The PCA was operating with two PCUs through- 
out the time range of interest here. However, only PCU2 was operating in the entire time span while 
PCU3 was active in the first 5.7 ks and PCU4 was operating in the remainder of the observation. 
Therefore, we accumulated spectral data collected with PCU2 only to avoid any calibration uncer- 
tainties. Given the fact that GX 301-2 is a bright X-ray source, we could obtain high signal to noise 
spectra even with one PCU. We accumulated phase integrated source and background spectra for 
the 12 time intervals shown in Figure [2] We modeled X-ray spectra in the 3-25 keV range, where 
the sensitivity of PCA instrument is the highest (Jahoda et al. 2006). 

We first fitted each spectrum individually with an absorbed power law plus a gaussian model. 
We find that this simple model fits all 12 spectra adequately well (xl ranging between 0.66 and 
1.23). We obtain that the hydrogen column density varies between 1.3x10^^ and 1.7x10^^' cm"^ but 
consistent with each other within errors. Similarly, the centroid energy and width of the Gaussian 
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Fig. 4. Power lax indices that are obtained by spectral modeling (filled circles) and the background 
subtracted light curve of the source in the 3-25 keV range (blue histograms). The horizontal bars 
indicate one spin cycle (i.e., 686 s) over which each individual spectrum is extracted. The vertical 
bars are Icr errors in power law indices. 



function (that accounts for Fe emission) do not vary across the accumulated spectra. Therefore, we 
performed a simultaneous fit of all 12 spectra by linking absorption, Gaussian centroid and width 
parameters but allowing power law index and normalizations of both parameters to vary for all 
spectra. The resulting best fit yields Nh = (1.6+0.2)xl0^^ cm"^ with a fit statistics of / degrees 
of freedom - 531.4/551. Note that the Nh we obtain here is consistent with what is determined 
for the orbital phase interval of 0.62-0.65 using BeppoSAX observations (La Barbera et al. 2005). 
We find that the spectrum of the source becomes softer starting with spin cycle number 4 and the 
softening trend peaks at the dip (see Figure |4]|. The spectral energy distribution is restored back to 
the pre-dip shape immediately after the dip and remains constant (within errors) in the spin cycles 
investigated here. 



4. Discussion 

We unveiled intriguing properties of the first observed X-ray dip from GX 301-2: (i) The X-ray 
pulsation from the source almost completely disappears for one spin cycle only during the dip and 
comes back to light in the spin cycle immediately after the dip, (ii) the phase-integrated spectrum 
of the source becomes softer right before and during the dip and it hardens again in the spin cycle 
immediately after the dip. Both properties point toward the fact that the mechanism for pulsations 
gradually turned off briefly, therefore, leading to the observed dip in GX 301-2. 

Although the observed dip in GX 301-2 is rather unique, other massive wind accreting X-ray 
binaries also exhibit a similar dipping behavior, in particular, VelaX-1 (Pspin = 283 s). In a long 
term hght curve taken wifli INTEGRL/lSGRl of Vela X- 1 , [Kreykenbohm eTall ( |2008] l observed 
several sudden drops in the count rate: similar to the dip discussed in this work, VelaX-1 was also 
pulsating normally at a flux level of ~250 mCrab, when the pulsations ceased in less than a minute. 
After about four spin cycles, pulsations resumed back to normal. During these dips in VelaX-1, the 
source was undetectable in 20keV to 40keV band by INTEGRAL /ISGRl. 
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Another such interesting dip was observed in a short RXTE observation of again VelaX-1 



(Kreykenbohm et al. 1999 1: at the beginning of the observation, VelaX-1 was not pulsating when 



the flux was at a very low level (37 counts/s/PCU in the 3-25 keV compared to 300 counts/s/PCU 
during normal bright state; see the inset of Figure |5]). Our re-analysis of these data (ObsID: 10141- 
01-03-00) shows that the source is detected in the PCA, but at a very low flux level and indeed no 
pulsations are present during this period. We extracted X-ray spectra from the time period when 
the source was in this dim non-pulsating state and from the bright pulsating state later in the same 
observation (see Figure|5]l. The spectral change between the two states is already evident when we 
rescale the spectral model of the bright pulsating state to the flux level of the dim phase spectrum: 
the spectrum of VelaX-1 is significantly softer during the dim non-pulsating phase. Therefore, 
GX 301-2 and VelaX-1 show the same behavior: they both exhibit short periods of time during 
which they are very dim, no pulsation is seen, the spectrum is significantly softer and recovery 
back to pulsing state is abrupt. 
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Fig. 5. Comparison of two X-ray spectra of Vela X-1 during the dip (filled circles, extracted from 
the time segment marked with dark gray in the inset) and in pulsing phase (empty circles, extracted 
time segment marked with light gray in the inset). The solid line is the rescaled model curve that 
describes the non-dip spectrum for comparison with the spectrum of the dip. 



Gradual reduction and disappearance of X-ray pulsation in GX 301-2 is suggestive of the fact 
that X-ray emitting regions of the source were blocked by the wind material ( [Kretschmar et al."| 
1999[ l. However, one would expect spectral hardening during the dip as softer X-ray photons would 
be absorbed. Since we observe significant spectral softening during the gap, it is unlikely that the 
dip is caused by stellar material moving into the line of sight. GX 301-2 is a source that has ex- 
hibited rapid spin up episodes ( |Koh et aT] |1997| l. As a result of an extremely rapid spin-up, the 
infalling material can be expelled via propeller effect. Cui (1997) reported the disappearence of X- 
ray pulsations in GX 1-1-4 and GRO J 1744-28 during low luminosity states and suggested that they 
are likely caused by the propeller effect. In case of the dip in GX 301-2, however, it is unlikely due 
to the propeller effect since the system would likely remain in the propeller regime much longer 
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than one spin cycle. Another possibility is that there could be a decrease and brief cessation of 
matter flowing onto the magnetic poles of the neutron star. In this scenario, the observed low level 
emission could be from the underlying polar region - a much larger area than the size of the accre- 
tion column which would be activated by recent inflow of matter and emission of which can easily 
be blanketed much brighter emission of accretion. It is still not straightforward to accommodate in 
this picture why the complete disappearence of the pulse lasts only one spin cycle. 
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